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Abstract 

The pion electromagnetic form factor with leading and next-to-leading twist cor- 
rections are studied in the framework of pQCD approach. We find that, at small 
momentum transfer regions, Sudakov effects cannot provide a strong enough sup- 
pression of the non-perturbative contributions coming from large transverse sep- 
arations, but at large momentum transfer region, non-perturbative contributions 
can be effectively suppressed. So pQCD approach can be applied at large momen- 
tum transfer region. At the energy region where experiment can access so far, i.e. 
Q < 4GeV, pQCD prediction can not be precise because there may be a quite 
large part of contributions coming from non-perturbative region. The calculation 
of power corrections may not be performed systematically in pQCD framework. 
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1 Introduction 



There is a general agreement that the standard approach of perturbative QCD (pQCD) 
is the correct theory in exclusive processes in the asymptotic limit Q ^ C)0 Although 
reasonable this theory seems in the ideal world, the test of this theory is performed only 
in the pre-asymptotic energy region. Up to now, the success of pQCD framework in 
exclusive process is very limited, such as in pion transition form factor where there is one 
hadron involved 0. The more hadrons involved in exclusive process, the less prediction 
power of the pQCD framework. It is found that the pion electromagnetic (EM) form 
factor contains substantial soft endpoint contributions which destroys the consistence of 
perturbative method 0. The mechanism of Sudakov suppression is introduced to suppress 
the soft endpoint contribution and a modified pQCD formula which includes Sudakov 
suppression is given in pion EM form factor It is claimed that pQCD calculation can 
still be self-consistent at about Q ~ 20Aqcd (2GeV for Aqcd = O.lGeV). 

During the past decade, there is no conceptual development of this modified pQCD 
approach. However the applicability of this method have met many theoretical prob- 
lems. Recently, this modified pQCD approach (or say pQCD approach for simplicity) was 
largely applied in exclusive B decays to treat the endpoint singularity . There has been 
much debate concerning the applicability of pQCD approach in B system 0. In 0, we 
investigated the reliability of pQCD approach m. B ^ -k form factors and find that the 
soft contribution coming from large transverse separations can be comparable with the 
hard contribution. This conclusion should be general for many exclusive processes, such 
as pion EM form factor, etc. . From the analysis in the light-cone sum rule approach, 
the soft endpoint contribution is found to be large at the experimentally accessible energy 
region [Q. So, the study of the reliability of pQCD approach in EM pion form factor is 
necessary and important. 

It is well-known that the leading order (LO) result for pion form factor with asymptotic 
leading-twist distribution amplitude (DA) is small, about | to | of the experimental 
data 0. The recent next-to-leading order (NLO) calculation is given in [^. Their 
conclusion is that NLO calculation is reliable only when Q > 5GeV. Other important 
corrections to the leading-twist LO result are higher twist corrections. The next-to-leading 
twist (twist-3 in our case) contribution is "chirally enhanced" power correction and it is 
most probably important. The study of twist-3 contribution to pion EM from factor was 
performed long ago [11]. The recent research using pQCD approach is given in |12|. Both 
of their results found a large twist-3 contribution even larger than leading twist result at 
intermediate energy region of Q < 5GeV. The large twist-3 contribution seems to destroy 
twist expansion. Since there is endpoint singularity in twist-3 contribution, one may also 
doubt the effectiveness of Sudakov suppression. 

In this paper, we will provide a systematic study of EM pion form factor in pQCD 
approach. Our main concerns are the reliability of applying pQCD method. Some new 
theoretical ingredients which are not considered in the previous literatures will be also 
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included in our analysis: intrinsic transverse momentum effects and threshold resumma- 
tion. We find that intrinsic transverse momentum in pion wave function are important. 
However, the large twist-3 correction shows that the power corrections may not be sys- 
tematically calculated in pQCD framework. 



2 Pion form factor in pQCD approach 

The pion EM form factor is defined by the following Lorentz decomposition of biquark 
current matrix element 

{n{P')\J,mn{P)) = {P + PXF^iQ') (1) 
where = J2 ^iQilfiQi is the electromagnetic current with quark flavor i and relevant 

i 

electronic charge Cj. The momentum transfer is = {P' — P)^ = ~Q^. We have 
restricted our discussion in space-like region. It is convenient to use light-cone variables 
in which P = (■^, 0, Oj_) and P' = (0, 6±). F^{Q'^) is the pion EM form factor which 
depends only on momentum transfer Q. The pion form factor at large momentum transfer 
Q provides information about the internal structure of pion. 

The basic idea of pQCD approach is that it takes into account transverse momentum 
and Sudakov suppression. The pion EM form factor is expressed as the convolution of 
wave functions V and hard scattering kernel Th by both the longitudinal momentum 
faction and the transverse impact parameter h: 

FTriQ"^) = J dx dy d% d%' V{x,b, jj,) V{y,b', jj,) TH{x,y,b,b',Q, jj,). (2) 

The wave function V{x, b, Q, /x) is given by: 

V{x, b, Q, /i) = exp[-S{x, b, Q, fi)]^%x, b) (3) 

where \E'°(a;, 6) is the soft part of pion wave function with \k_^_\ < 1/b 

^°(x,6) = 0(x,l/6) + O(a,(l/6)) (4) 

The above equation is valid for small b. When Sudakov suppression is strong, there is 
only small b contribution and the approximation of 6) by distribution amplitude 

1/6) is valid. But at a few GeV region, this approximation is questionable. 
The factor exp(-S) in Eq.(|^) includes the Sudakov logarithmic corrections and renor- 
malization group evolution effects of both wave function and hard kernel, 

S{x, b, Q, /i) = s{x, b, Q) + .(1 - X, 6, Q) - lln^^^;^^ (5) 

where /3i = with = 3. The Sudakov exponent 6, Q) is calculated up to next- 

to-leading-log (NLL) accuracy. Its explicit formula can be found in fl^. The exponent 
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s{x,b,Q) is obtained under the condition that xQ/\^ > 1/6. For small b, there is no 
suppression, so s{x, b, Q) is set to zero for xQ/a/2 < 1/6. 

The study of distribution amplitudes beyond leading twist is expanded in the con- 



formal spin [jT4| . The light-cone distribution amplitudes of pion are defined in terms of 
bilocal operator matrix element 
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(6) 



a/3 



where /^r is the decay constant of pion and x is the longitudinal momentum fraction of 
quark in pion. The parameter /i^ = m^/irriu + md) for charged pion. 0,^, 0p and (pa are 
the twist-2 and twist-3 distribution amplitudes, respectively. The twist-3 terms contribute 
power corrections. At the experimental accessible energy region, the chirally enhanced 
parameter = /Xtt/Q ~ ^'(l) is not small. 
The momentum projection for pion is [0]: 



'75 - yUp75 (pp - ify^iu 



p ■ p 6 



d 



6 dk±iy 



(7) 



0/3 



where 0^ = 



The final formula for pion EM form factor in pQCD approach is: 
FAQ') = Y^/'Q' I' d^dy bdb b'd' X [y-VM'PM + 4 



+ (l+y)P,[x]^ + 2Pp[x]^'^f^] 



6 



6 



X H{x,y,b, b',Q,ii) 



y'Pp[x]Vp[y] 



p{2) _^ p{3)_ 



where F^^^ , F^^^ represent twist-2 and twist-3 contributions in pQCD approach respec- 
tively, and H is given by 



H{x, y, b, b', Q, ji) = Ko{y/xy Qb) 



e{b-b')Ko{VVQb)Io{VVQb') 



+eib'-b)KoiVVQb')Io{VVQb) 



where [x] = {x,b,Q,fi), [y] = {y,b',Q,fi). The wave functions of V-^.Vp and Va can be 
obtained from the relevant wave functions ^t,, \E'p and through Eq.(|) and Eq.@. Kq 
and /q are the modified Bessel functions. The choice of renormalization scale parameter 
H is taken as the largest momentum scale associated with the exchanged virtual gluon in 
the longitudinal and transverse degrees. 



/i = max(v/xy(5, 1/6, 1/6') 



(9) 
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The above choice avoids the Landau pole in couphng constant ^^(yu) at /i = Aqcd if x 
and y are smalL 

In [0 , the transverse momentum associated with virtual fermion lines is neglected and 
the hard kernel involves only a single impact parameter h. In our formula of Eq. (P) , there 
are two parameters h and h' . The errors caused by the neglect of transverse momentum 
in the propagator of fermion lines is small if Q > 4GeV, but it can be be about 20% 
when Q = 2GeV. So, we retain the transverse momentum in the fermion lines without 
approximation. The transverse momentum k'j_ in the numerator are neglected because it 
is power suppressed compared to Q^. We checked this assumption and found its effects 
are really small in our case. 

In the asymptotic limit, the distribution amplitudes (pn{x) = 6xx, (f)p{x) = 1 and 
4>cr{x) = Qxx. Neglecting the transverse momentum in both hard kernel and wave func- 
tions, Eq.(RD becomes 
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1 dx 



X 



1 dx 



X 



(ppix) 



Ff+Ff) (10) 



where Fj^^ represent twist-2 and twist-3 contributions to Ft^^Q"^) in the standard 
approach respectively. The twist-3 contributes to correction and it is power sup- 
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pressed by a factor ^ compared with the leading twist contribution. Because 0p(x) is a 
constant at endpoint, twist-3 contribution is logarithmically divergent at x = 0. We use 
the effective scale /le// in Eq. ([ToD because the NLO calculation depends crucially on the 
choice of /x. 



3 Numerical results and discussions 

There are only two parameters Aqcd and in the hard kernel. They are chosen as 
Aqcd = 0.2GeV and ix-,, = 2.0GeV. The distribution amplitudes are taken as their 
asymptotic limit form. We do not use C-Z distribution amplitude for discussion since 
this model of distribution amphtude are concentrated at the endpoint where perturbative 
analysis is not reliable. 

At first, we discuss a problem which is not investigated after the paper of the 
consistence between the resummed formula and the standard formula. The modified 
pQCD approach uses a b-space resummation formalism in which the Sudakov double 
logarithms are resumed to all orders. There is a difficulty of matching the resumed 
formalism and the fixed order predictions in the standard approach. From the theoretical 
point of view, the modified pQCD approach should be consistent with standard approach 
when Q is large enough. Because of the existence of endpoint singularity at twist-3 level, 
twist-3 contribution cannot be calculated in the standard approach. We can compare the 
predictions of pQCD approach and standard approach in leading twist, i.e., F^^^ and F^^^ 
to check the consistence. 
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The physical quantity F„ does not depend on the choice of the scale parameter n if the 
calculation can be performed up to infinite orders. However, in practice the calculation 
can only be made perturbatively at finite orders. To make the perturbative expansion 
reliable, the scale parameter should be chosen in such a way that make the higher order 
corrections as small as possible. In pQCD approach, fi = maxl^/xyQ , 1/6, 1/b'). While in 
the standard approach, it is a free parameter. We use the different choices proposed in 
literatures: fieff = Q,XQ,X^/^Q,e-^^^XQ where X =< x >= 1/2. Fig.l shows the Ff) 
in the standard approach with different choice of /i and Fj^^\ One can see that the best 
choice of /i in the standard approach is /i = Q/2 or X^/^Q. By using this choice, pQCD 
approach is consistent with the standard approach when Q > lOGeV. For the momentum 
transfer Q < lOGeV, the prediction of pQCD approach is generally smaller than that in 
the standard approach. This difference is caused by the transverse momentum effects and 
Sudakov suppression. In [|T^, it is found that the choices of /i = XQ, X^^'^Q reduce the 
NLO corrections significantly and reliable predictions can be obtained at lower values of 
Q. Both these results show the importance of choosing the scale by the interior dynamics 
of the process. 




0.05 ^ : r I - 

2.5 5 7.5 10 12.5 15 17.5 20 



Q(GeV) 

Figure 1: The comparison of twist-2 pion form factor in standard approach and pQCD 
approach. The solid line is for pQCD approach, while the dashed lines for standard 
approach: (a), (b), (c) and (d) for fx^ff = Q, XQ, X^/'^Q and e~^^^XQ, respectively, 
where X =< x >= 1/2. 

Now we discuss the reliability of pQCD approach. The basic idea of pQCD approach 
is to use Sudakov form factor to suppress the long-distance contributions coming from 
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large transverse separations. A reliable perturbative computation should satisfy that 
most of the results are coming from small impact parameter b. In order to study the 
impact parameter b dependence of form factor F^, we introduce a cut-off b"^ in the impact 
parameter space in the integrals of Eq. (P) by Jq db. Similarly, the impact parameter b' 
dependence can also be performed. It is convenient to use Q^F as the physical quantity 
for discussion because Q^F is nearly a constant for large Q. The results are plotted in 
Fig. 2. We can see that the results are saturated before 6^ approaches to 1/Aqcd, which 
means large separation is suppressed by Sudakov effect. The saturation point is closer 
to 1/Aqcd for smaller momentum transfer. For Q = lOGeV the saturation point is at 
about 1.5GeV~^, which is far from the end point 1/Aqcd- This means that almost all the 
contributions come from the short- distance region. For Q = 6GeV, the saturation point 
is at about 2.5GeV~^, which shows that some non-perturbative contribution emerges, 
but it is still not large. But for Q = 4GeV and 2GeV, the saturation points are at 
3.5GeV^^ and 4.0GeV~^, which are quite close to the end point 1/Aqcd- There are 
substantial contributions coming from large transverse separations b > 0.5/Aqcd for Q < 
4GeV. Sudakov suppression becomes weak for small Q, and non-perturbative contribution 
becomes large. To show directly how the non-perturbative contribution becomes large at 
small momentum transfer, we show the dependence of Q'^F{Q'^) in Fig. 3. We see 
that there are quite large contributions coming from the region Q;s(/i) > 0.5 for smaller 
values of Q"^: 34% at = 4GeV^ 22% at = lOGeV^ 18% at = leGeV^. For 
> 25GeV^ the contribution for ^^(/i) > 0.5 becomes smaller than 10%. 




Figure 2: Dependence of Q^F on the cut-off b^ aX Q = 2GeV , AGeV, 6GeV and Q = 
lOGeV. 

Is Sudakov suppression effective at experimental accessible energy region [Q = 1 — 
4GeV)? From the above results, it seems that Sudakov suppression is not strong enough 
for Q < 4GeV. Sudakov form factor Q-si^'>')-s{x,b) ^^^y suppress the region of 6 ~ 1/Aqcd 
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at a few GeV region. In the derivation of b-space resumed Sudakov form factor, it 
requires the condition Q ^ 1/6 ^ Aqcd- The condition 1/6 ^ Aqcd is required because 
a rehab le as(l/6) expansion needs smaU h. This is possible only if Q is large enough and 
Sudakov form factor suppress all the large h contributions. The derivation of Sudakov form 
factor is a difficult problem in QCD because perturbative expansion is not meaningful at 
long-distance. This is unlike the case in QED. Thus, Sudakov form factor depends on 
the infrared cut-off. Thus, b-space resummed Sudakov form factor is not the unique 
choice. Moreover, the extrapolation of asymptotic form into the small momentum region 
is not under well control in QCD. The predictive power of perturbation theory will be 
decreased largely by this extension. The model-dependent nonperturbative effects have 
to be included in order to extrapolate the perturbative analysis to the large h region with 
the cost that the predictive power of perturbative method is decreased. In the pion EM 
form factor, the non-perurbative contribution at small Q can be at the order of 30%. This 
means that Sudakov suppression is not strong enough for small momentum transfer. 
Other mechanisms, such as intrinsic transverse momentum effects IqII and threshold 



resummation can suppress large h contribution and endpoint contribution. We in- 
vestigate whether these effects can provide a reliable perturbative analysis. About the 
intrinsic transverse momentum effects, it is nonperturbative. The estimate of this effects 
must be model-dependent at present. The suppression of the intrinsic transverse momen- 
tum effects is larger than Sudakov effects at the order of a few GeV. Only for very large 
Q, Sudakov suppression becomes strong and the intrinsic transverse momentum effects 
can be neglected. We incorporate these effects to suppress the nonperturbative contri- 
butions. The formula of including these effects can be found in [0. Here, we do not 
present them for simplicity. Fig.4 shows the numerical result after incorporating these 
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effects. The result is seriously suppressed by intrinsic transverse momentum effect in pi- 
ous' wave functions, and slightly suppressed by threshold resummation, which shows that 
intrinsic transverse momentum effect is important in pion EM form factor. We should 
include them in our analysis. Fig. 5 shows that the perturbative behavior is improved 
largely. The saturation points in 6-space moves to the 6^0 direction apparently for 
both small and large value of Q. The contribution of as{ij) > 0.5 is largely reduced. 
For > 20GeV^, the contribution of as{ij) > 0.5 is suppressed below a few percent. 
However, for ~ 4 — lOGeV^ there arc 10% to 20% of the contributions coming from 
the region ^^(/i) > 0.5. Non-pcrturbative contribution for small Q region is still non- 
negligible. So pQCD can not predict pion EM form factor precisely in small momentum 
transfer region. For large momentum transfer region, the non-perturbative contribution 
is small and pQCD can be applied. 



. 6 p. 
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Q2(GeV2) 

Figure 4: Dependence of Q^F on Q^. The solid line is the result including intrinsic 
k±_ in pion and threshold resummation effect; the dotted line is without intrinsic k±_ and 
threshold resummation; the dashed line is with intrinsic k±_ effect but without threshold 
resummation. 

^From the phenomenological point of view, we can use pQCD approach to estimate the 
hard contribution. The pQCD approach is valid in perturbative regions where the strong 
coupling constant as{^J) be small. In order to estimate the hard contribution, one should 
make a criterion that a perturbative contribution should satisfy although it is very hard 
to define such a criterion quantitatively. In general perturbative contribution should come 
from the region with small strong coupling constant a^. To make our analysis proceed 
numerically, we set a criterion for perturbative contribution: < 0.5. This criterion can 
not be understood as an absolute one, it is only indicative. A more stronger criteria may 
be more reliable but may lose some hard contributions. 

Table 1 gives the numerical results which satisfy ag < 0.5. From it, it is seen that 
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Figure 5: Dependence of Q'^F on: (a) the cut-ofT 6"^ and (b) Q^, with intrinsic transverse 
momentum effects and threshold resummation effects included. 

Table 1: The dependence of the hard contribution of pion EM form factor 
Q'^F on Q. The rows "twist-2" and "twist-3" represent twist-2 and twist-3 
contributions. The "total" represents the sum of twist-2 and twist-3 contribu- 
tions 



g2(GeV^) 


1 


2 


4 


6 


8 


10 


16 


25 


twist-2 


0.024 


0.042 


0.059 


0.066 


0.071 


0.073 


0.077 


0.078 


twist-3 


0.36 


0.30 


0.22 


0.17 


0.14 


0.12 


0.086 


0.060 


total 


0.38 


0.34 


0.28 


0.24 


0.21 


0.19 


0.16 


0.14 



the perturbative contribution is about 0.16 for Q > 4GeV with LO twist-2 and twist-3 
corrections. Compared with the experimental data for IGeV^ < < 6GcV^, the results 
in Table 1 arc not too bad. Considering the large error bars in the present experimental 
data, we cannot conclude that soft contribution in pion EM form factor at small region 
can be neglected. About 20% percent soft contribution at small momentum transfer region 
is still allowed. For Q"^ — IGeV^, the soft contributions is larger than 30%, thus pQCD 
can not be applied in this energy region. On the other hand, the results in Table 1 can not 
be viewed as the self-consistent prediction of pQCD because we set as < 0.5 in getting 
them. Only for large Q^, the results with the constraint ag < 0.5 can be very near to the 
full pQCD prediction. As becomes large soft contribution decreases very fast, hard 
contribution seriously dominates. Table 1 also shows that contribution of twist-3 is large 
for small Q^, it becomes to be smaller as being larger. 

About power correction, it may contain both hard and soft contributions from a gen- 
eral point of view. In pion form factor, we find a large twist-3 contribution which is 
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enhanced by the endpoint region in pQCD approach. The large twist-3 contribution 
means that twist expansion may be questionable. Because Sudakov suppression is weak 
at low momentum transfer Q, the soft contribution can not be suppressed effectively. 
Whether twist-3 contribution is large or small should be checked by other complementary 
non-perturbative method in order to guarantee the validity of twist expansion. 

At last, we briefly compare our perturbative analysis with the prediction from the 
QCD sum rule. There is a major problem about OPE (operator product expansion) in 
the early three-point QCD sum rule so that QCD sum rule is hold for Q < M where M 
is the Borel parameter which is at the order of 1 — 2GeV^ §]. At this energy region, the 
factorization of the hard and soft dynamics is generally difficult or impossible. In the light- 
cone sum rules, the expansion is carried out by twist of hadrons so that comparison of the 
perturbative framework and light-cone sum rule becomes possible. A systematic analysis 
of pion form factor including the radiative corrections in the light-cone sum rule approach 
is provided in [§]. They found that the soft endpiont contribution is dominant and the 
hard contribution is very small. Our analysis also shows LO twist-2 hard contribution 
is very small. But the calculation of the hard part of twist-2 contribution is different 
in these two approaches in principle. In pQCD approach, the momentum of exchanged 
gluon is determined by the longitudinal and transverse momentum of the quarks in the 
two pions. The LO calculation is at the tree level. In the light-cone sum rule, the LO 
order is soft endpoint contribution and the NLO radiative corrections contains hard and 
soft contributions. There is a certain seeming correspondence between the single-gluon 
exchange diagram in pQCD approach and a diagram of radiative corrections in the light- 
cone sum rule approach. However, they are not equivalent. In general, the loop corrections 
are basically different from tree level calculation. In loop corrections, the momentum is 
arbitrary and the result has ultraviolet divergence and infrared divergence. 

In conclusion, in large momentum transfer region, say > 20GeV^, non-perturbative 
contribution is small, EM pion form factor can be self-consistently calculated by pQCD 
approach. While, in small momentum transfer region, < lOGeV^, Sudakov suppression 
becomes weak, non-perturbative contribution can not be completely suppressed. There is 
always at least 20% soft contribution left in pQCD prediction if we use the soft criterion 
as > 0.5. This soft contribution breaks the self-consistence of pQCD approach. The 
solution of this problem should be beyond the perturbative framework. The reliable 
calculation of power correction requires non-perturbative methods. 
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